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The extended Skyrme-Faddeev model possesses vortex solutions in a (3+1) dimensional Minkowski space- 
time with target space CP^. They have finite energy per unit of length and contain waves propagating along 
vortices with the speed of light. We introduce various types of the potentials which correspond with holomorphic 
solutions of the integrable sector and also with several numerical solutions outside of this sector. The presented 
solutions constitute a strong indication that the current model contains large class of solutions with much wider 
range of coupling constants than the previously known exact solution. 
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I. INTRODUCTION 


II. THE FORMULATION OF THE MODEL 


The Skyrme-Faddeev model is an example of a field theory that 
supports the finite-energy knotted solitons. The significance of this 
model has increased noticeably when it has been conjectured that the 
model can be seen as a low-energy effective classical model of the 
underlying Yang-Mills theory Ql]. Similarly to many other models 
13] the classical soliton solutions of the Skyrme-Faddeev model can 
play a role of adequate normal models useful in description of the 
strong coupling sector of the Yang-Mills theory. The exact soliton 
(vortex) solution of the model has been found within the integrable 
sector d. Such a sector exists in the version of the model that is 
an extension of the standard Skyrme-Faddeev model obtained by in¬ 
cluding some quartic term different to the Skyrme term. The study of 
the extended models have been originally motivated by the results of 
the analysis of the Wilsonian action of the SU (2) Yang-Mills theory 
d- It has been shown that also in the case of the complex projective 
target space CP^ the extended Skyrme-Faddeev model (in which it 
has been imposed some special constraints for the parameters of the 
model) possesses an exact soliton solutions in the integrable sector 
SIS . It has not been clear until now if the presence of the solutions 
in the considered model is related to the particular choice of the cou¬ 
pling constants like in the case of the exact solution or it is rather a 
general property of the model. In order to answer this question one 
needs to construct some other solutions than the exact ones. The aim 
of this paper is to investigate the existence of the solutions of the 
model inside/outside the integrable sector, especially in absence of 
some particular relations between coupling constants. We show that 
such solutions exist. The key role is playing by the potential which 
usually work as a stabilizer for the solution. In the present paper the 
potential appears in the context of exact holomorphic solutions, as it 
was presented in the case of CF' related model aa, and also for 
solutions from a non-holomorphic sector. We conclude that the exact 
solution appears as a particular solution belonging to the wider class 
of solutions of the model. 

The study of such models is promising and could be important 
for understanding some aspects of the strong coupling sector of the 
Yang-Mills theory. 
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The Skyrme-Faddeev model and its extensions on the CF* target 
space are usually expressed in terms of the real unit vector n. The 
dimension of a target space is simply related to the number of degrees 
of freedom of the model. For instance, the model with the CF* target 
space has only two independent degrees of freedom. In order to add 
more degrees of freedom one can consider some higher dimensional 
target spaces. The target space (coset space) in the case of some 
higher dimensional SU (N) Lie groups, i.e. N > 2, can be chosen in 
several nonequivalent ways. 

Recently it has been proposed some formulation of the extended 
Skyrme-Faddeev model on the CP^ target space (1 . The coset space 
CP^ = SU{N + l)/SU{N) 0 (7(1) is an example of a symmetric 
space and it can be naturally parameterized in terms of so called prin¬ 
cipal variable X{g) = g(T(g)^*, with g e SU (A+ 1) and a being the 
order two automorphism under which the subgroup SU (N) 01/ (1) is 
invariant i.e. o(h) = hfoi hiz SU(N) 0 (7(1). The principal coordi¬ 
nate X (g) defined above satisfies X (gh) =X(g). 

We shall consider fhe field fheory in (3 H- 1) dimensions defined 
by the Lagrangian 
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where M is a coupling constant with dimension of mass whereas 
the coupling constants j3, y are dimensionless. The first term 
is quadratic in Y and corresponds with the Lagrangian of the CP^ 
model. The quartic term proportional to e^^ is the Skyrme term 
whereas other quartic terms constitute the extension of the standard 
Skyrme-Faddeev model. The novelty of the model 0 comparing 
with that introduced in @1 is the presence of the potential V. In re¬ 
cent studies, several potentials have been introduced for the planar 
Skyrme-type model (^- [TTIl . It was shown that the extended CF* 
Skyrme-Faddeev model in (3+1) dimensions possesses some non- 
holomorphic solutions that do not belong to the integrable sector (3). 
Since the extended Skyrme-Faddeev model on the CP^ target space 
also possesses the integrable sector as well as the exact vortex solu¬ 
tions the natural question is if there exist any solutions that do not 
belong to the integrable sector? In similarity to the paper j3] we 
study such a possibility in the presence of the potential. As it has 
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been explained below such solutions can be obtained numerically for 
some choice of the potential. 


A. The parametrization 


Let us shortly discuss the parameterization of the model. Ac¬ 
cording to the previous paper |3] one can parametrize the model in 
terms of N complex fields m,-, where i= 1,... ,N. Assuming (A+ 1)- 
dimensional defining representation where the SU{N+ 1) valued el¬ 
ement g is of the form 
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and where A is the hermitian N x A-matrix 

UiU*i 




which satisfies A-u = u and ■A = u' 


l + tJ 

The principal variable A (g) = gCT(g)^* takes the form 
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and the Lagrangian Q reads 
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where the symbols C^v and are defined as follows 
Cjiv '■= M^T]iiv ~ ^ 

+{ye^ - l)-r^v + + 2 )Tv;j j, 

• d^U~ [dyu' ■ u){u^ ■ dfi 


'^IIV — 




(4) 


(5) 


A variation with respect to u* leads to the equations which can be 
cast in the form 

• u)d^{Cjiyd'^Ui) 
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where we have already multiplied the resultining equations by in¬ 
verse of Afj i.e. AJj^ = „ i8ij + UiU*j). We shall discuss some 

examples of the potential in the further part of the paper. In the sim¬ 
plest case when the potential is a function of absolute values of the 
fields V(|ni|^,...,|n^p) the contribution from the potential becomes 


N r 
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We introduce the dimensionless coordinates (r,p. (p,z) defined as 

x^ = rQt, jc'= ropcostp, = ropsintp = rgz (7) 

where the length scale rg is defined in terms of coupling constants 
and i.e. 

° M2e2 


and the light speed is c = 1 in the natural units. The linear element 
ds2 reads 

^ rQ{dt^ - dz^ - dp^ - p^dqp-). 

The family of exact vortex solutions has been found for the model 
without potential = 0 where in addition the coupling constants 
satisfy the condition /3e^ + ye^ = 2. The exact solutions have the 
form of vortices which depend on some specific combination of the 
coordinates i.e. one light-cone coordinate x’ +x® and one complex 
coordinate x* + ix^. The functions n,(x^ +x®,x* -I- ix^) satisfy the 

zero curvature condition d^Uid^Uj = 0 for all i,j = and 

therefore one can construct the infinite set of conserved currents. 

We shall consider the following ansatz 

uj = ( 8 ) 

where V'(w') is a real function of the light-cone coordinate and fj{p) 
are real-valued functions. The constants n; form the set of integer 
numbers and kj are some real constants. The holomorphic solutions, 
which belong to the integrable sector, is of the form 

Mp) = Cip"' (9) 

where c, are some real (in general complex) free constants. We define 
two diagonal matrices 

A =diag(ni,...,nAr), (T = diag(ki,... .fc/v). (10) 


in order to simplify the form of some formulas below. In matrix 
notation the ansatz reads u = f{p)exp[i{?.(p + (7\lf(w))] where w is 
either z + t oi z — t. The expressions have the following form 


tpp ^ e(p) = - ^ (/^•/)(/^./')] 

= (0{p) = 

T(j)p = —T^pip = i?(p) 

^(p)=(/^.A./)(/^./)] 

T,p = -Tpr = {id,\l/)^{p), X,p = -Tp^ = (id-,^f)^{p) 

^(P) = [tjV'^.CT./- (/.(T./)(/^./)] 

= {d,\i/)ri{p), = (d,}j/)rj(p) 

P(P) = 

Xtt = {d,\l/fx{p), ^u. = (dzwfxip) 

= % = {d,\lf){dz\l/)x{p) 

Xip) = 

where derivative with respect to p is denoted by — and T stands 
for matrix transposition. The equations of motion written in dimen¬ 
sionless coordinates take the form 


( 1 +/^./) 
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for each k= \ where we have introduced the symbols Cfjv '■= 

and also := The components Cuv which appear 

in the equations of motion read 

Cpp = —1 + (/3e^ — 1) ^6 + + (2ye^ + 1)6 

Cff = + P^ipe^ — 1) ^6 + T + 1)® 

c^p =-Cp^ =-3i^. (12) 


B. The energy 


The Hamiltonian density being a Legendre transform of the La- 
grangian density 0 is defined as follows 


:= — doUi+ dpfU*-^. (13) 

O OQUi O OQUj 

The resulting Lagrangian and Hamiltonian densities taken for the so¬ 
lution 0 depend only on the coordinates p and either z + t oi z — t. 
For the Lagrangian density one gets 
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where a term proportional to is just the CP^ Lagrangian and 
terms proportional to ye^ — 1 and ye^ + 2, vanish for the holomorphic 
solutions 0 since the constraint dpuyuj = 0 leads to the relation 

f'jiP) = jfjiP) resulting in equalities e = -^ = ^ondy^. 

The Hamiltonian density GS considered for 0 contains several 
terms. We present it in dimensionless form (in the unit of —M^e^/4) 


7 

(15) 

7=1 

where the components M'j are given by 

,*i—( 9 +“). M.-2yyz. 
jti-y^v^-2)(e+y 

We have split the Hamiltonian density in order to make explicit the 
terms that were present in the earlier study of the holomorphic vor¬ 
tex solutions i.e and The term Ma, and also the po¬ 

tential term Mq were absent in previous considerations. They did 
not appear due to the constraint imposed on the coupling constants 


Pe^ + ye^ = 2 and also due to the integrability condition. The terms 
are zero for the holomorphic solutions. Note that for static 
(w-independent) vortex solutions, reduce to zero and 

then, only the terms could be meaningful for the Der¬ 

rick’s scaling argument. 


C. The topological property 

According to the discussions in in and also in @], we can define 
the topological charge in the present model. The field u; provide a 
mapping from plane into CP^. However, for the finiteness of 
the energy, the field goes to a constant. Then the plane should be 
compactified into and the solutions define the mapping CP^ 

which is classified into the homotopy classies of K 2 {CP^). There ex¬ 
ists a theorem describing in 111 , MG/H) = 7ti {H)g where 7ti {H)g 
is the subset of 7ti (H) formed by closed paths in PI which can be 
contracted to a point in G. Thus, in the present case, the topological 
charges are given by 

;ri(5(/(iV)®(/(l))sj;(^+i). (17) 

As discussed in 111,11, the topological charges are equal to the 
number of poles of u,-, including those at infinity. And then, it can 
be obtained as 


Qtop — Umax T |77minl (18) 

where the highest positive integer in the set n;,i = 1,2,••• ,N and 
Umin is ths lowest negative integer in the same set. 


III. REDUCTION TO THE INTEGRABLE CP^ SECTOR 


In H, one of us have claimed that for some special choice of 
the potential V there exist analytical solutions of the CF* Skyrme- 
Faddeev type model for all topological charges. It turns out that for 
a CP^ version of the extended Skyrme-Faddeev model there exists 
sectors such that the model reduces to the CP* version. One can 
expect that when reduction occurs the model possesses holomorphic 
solutions for the appropriate choice of the potential. In current sec¬ 
tion we shall study this problem in details. 

We are interested in a case such that the quartic term proportional 
to Pe^ + y — 2 does contribute to equations of motion and solutions 
of those equations are of the form 

uk = cjp«te'(«i'P+<^V'(H')) (19) 


with Cl; being some real constant parameters. The form of solution 
solves the zero-curvature condition dpuid^uj = 0. The problem we 
have to face is in fact an inverse problem i.e. we shall derive the form 
of the potential for a given solution. Let us observe that the equations 
of motion for the solution G1 simplifies a lot taking the form 


-2 ^(p) ^ p2 
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- £ 
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5V , ^ . .2 SV 
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where the coefficients are functions of integers n; 

A; =«?(«;- 1 ) 

Bij = nf{ni-2nj - 1) + ^{rii - nj)^{ni+ nj - 1 ) 

1 2 

Cijk = {ni + nj-lrik-l). 

One can get rid of the denominator (1 + • n)"^ on both sides of the 

field equations substituting 


SP 


V = - 


w 


(1 + nf • ti)4 

The resulting set of equations i= I,... ,N has the following form 


/12 p4 


R{p) = 


dW . ^ , 2 


5|t(i|2 ' S\uk\ 


41V. (21) 


The function IT must have such a form that the complete system 
of N equations (ED holds. Such a function IT, satisfying the set 
of equation with arbitrary integers { 111 , 112 , ■■■ pin), would be a true 
generalization of the problem of exact solutions to the CP^ case. 


A. The potential 

Instead of solving this (still open) question we shall study the case 
of reduction with only one non-zero integer n. We assume that the 
first K integers n; are equal nj = 112 = ■ ■ .tif; = n and rest of them 
vanish n/c+i = ... = «^ = 0. The set of reduced equations ED takes 
the form 

I = » , 22 , 

\ s&+Ef.,i..r-sgp-« = 0 

where the first subset of d 22 b is labeled by i = 1 ,..., A" and the second 
subset by i = A + 1,... ,N. The solutions m; for c\ = ... = ck = c 
satisfy 


where the parameter a is a free constant. Considering that 

P and Vf , In; 12 = KaP one can show that the second sul 

\uj\- ^J=l' j' 5|«j|- 

set i = A-f 1 ,..., A of equations l l 22 b reduces to the following one 


SQ 


8Q 




+ (Aci:-4)e = 0. 


(25) 


The equations l l25b considered for |m/| 2 = c? hold for any i = KP 
1,..., A if the function Q is such that 


1 SQ 




= 8 


(26) 


where 5 is another free constant. Consequently the equations l l25t 
reduce to the relation between constants a and S 


5 = 


4-Ka 


I + IIa 


^2 • 


(27) 


~j=K+l^j 

For the first subset i= 1,.. .,K of the equations ll22b one gets 


y ^ \ 


Ct ^r-1 

j^+Ka+ Y, 


j=K+\ 


jl SQ 
Q S\uf 


PQ 


where the rhs of this formula is taken for \ uj\^ given by l l23b . It follws 
that the last formula in fact became 




|m |2 


PQ- 


where P = |n|2^“, Q = 2(4... ,4). In the last step we have 
made use of the relation (1271) . Putting all results together we obtain 


|m,|2 =: |m| 2 = c2p2« for i=\,...K 
\ui\^ = cf for i = K+l,...N 


(23) 


One can conclude from the last equation that both free constants a 
and S have to be fixed by 6 : = cc and 5 = 5. In such a case one gets 


which implies that the function R{p) in ( 122b simplifies to the form 

PiP) = ^K^ (i+ E [a|n|2-5|n|4l 
^ V k=K+l J 


- 5(|n|2)H = ^ [a|n|2 - 5|n|4l 

and all equations take the form of relation between constants which 
fixes value of p 2 in terms of other constants 


where the coefficients a, 8 read 


a := 


2 

K 



5 := 


2 


i + L 


N 

k=K+l 



(24) 


The problem can be solved by the method of separation of variables. 
For this reason we consider the function IT in the form of a product 
of two functions P and Q 

W = P{\ui\^,...,\uK\^)Q{\uK+\\^,...,\ur^\'^) 


where the function P is a product 


= 32{l5e^ + Ye^ -2)c«n^K^ 


• , 4 ) 


Q{c 


(28) 


/£■+!’ 


The only condition that the function Q has to satisfy is that given 
by ( 126b with 5 = 8. An example of the function Q is given in next 
subsection. 

We conclude from this section that for the potential 


{Uf=iW,\^rQ{WK+il---,\m\^) 


P{M\...,\uK\^)=fl\u,\^" 


;=i 


where Q satisfy the condition ( 126b in the case of reduction to CP* the 
model possesses holomorphic solutions in the sector jSe^ + ye^ 7 ^ 2 . 
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An example of the exact potential 


B. The energy of an exact vortex configuration 


Let us consider the function Q in the form 


Q = 


n (l"d ) "‘S(l“tr+ih---,l“ivr) 

l=K+l 


In such a case the condition l l26t i.e. 


= 8 became 


d\uj\^ 


Ini' 


nc^ 2 


We shall assume the function S in the form of linear combination of 
|i<;P i.e. S = £/ — which reduce the last condition to 

the set of equations 


N 

^ E 

l=K+l 


+n I 1 + E C; 
l=K+l 


1+ E c?+('^+i) 

l=K+l 


CjSSj = 0 


that have to hold for any j = K . ,N. All coefficients propor¬ 

tional to c? as well as free terms must vanish independently. Finally it 
leads to two equations what suggest parametrization containing only 
two variables a and b. Taking c^jfSj = ac? -f we reduce the set of 
equations to the following one 

(LiL/r+icf ^n)fl+(n+l)(iV-Ai)fi= (n+1)^/ 

(it/r+l cli)a+{N-K-n)b = s^ 

that have solutions in the form 

( „ - _ ("+ 1 K 


V -(n+l)+(A'-A'+l)(l+L/Ls:+iC^) 

For a particular choice of i.e. 


'=-(«+i)+(iv-ff+i)(i+ X ch 


(29) 


l=K+\ 


solutions reduce to a = n + 1 and ^ = 1 + L/La'+i follows that 

1 


— ^ + 1 + —2 I 1 


N 

- E 

l=K+\ 


(30) 


The set of coefficients ,..., determines the constants and 
iMj. From the physical point of view the inverse problem is more 
interesting i.e. when the potential parameters are free constants. In 
such a case one has to invert the relations between cj and fSj. The 
parameter .e/ is not independent constant since it is determined by 
values of constants Cj. The function Q gets the form 


Q = 

which implies the formula 


n 

l=K+\ 


l=K+\ 


(31) 


q{ck+\t■ ■ ^cn)—+ E n ‘'i 

j=K+l l=K+l 

Plugging this result to l l28b one gets the relation 


In this subsection we present exact expressions for the energy of 
the vortex configurations being analytical solutions in the model con¬ 
taining potential. To make formulas less complex we shall consider 
a simplified case such that there is only one free constant Ck = c for 
all /: = and q = 1 for A; = AT-f 1,... ,N. It follows that 


2 / 1 

«:= -p 1 -- 

K \ n 




(32) 


The energy per unit of length of the vortex is a sum where 

contributions are defined as 


Et = 27 iM^ [ dppM- 

Jo 


The first term £1 is purely topological and therefore it is proportional 
to \n\ 

El = 9,nhf\n\. 

The energy £2 can be cast in the form of the sum 


£2 = 8 ;rMV^ 


1 


N-K+l 


K 


\n\c« L 

where the coefficients aj depend on parameters ki 


Y,ajlj(n) 
-> f=l 


K 


K K 


«1 := 


1^2 


vA' pi yK pi yK yN 

^ . Li=K+l'^i , ^Li=lLj=K+l^''^J 

02 ■= —— . H-77- 2 - 


^3 • = 


N~K+l ' K K{N-K+\) 

yN pi yN yN e.e. 
Li=K+l'^i U=K+\i~j=K+l'^i'^J 


N~K+l (N-K+lf 
and expressions Ij{n) stands for some integrals. One gets 

I 

y- 


/,(«) := r 

Jo 


dy 


;r(«+l) 


(i+y )2 


n = - 1 ,- 2 ,- •• . 


The integral /] (n) diverges for positive values of n. The integral I 2 (n) 
reads 




+y)^ 


n = ±2,±3,--- . 


where divergence occurs only for n = ±1. The last integral is of the 
form 




'!• ^ 7t{n — 1 ) 

+ y )2 nsin(f) 


n = 1 , 2 ,--- . 


One can see that there are not values of n such that all integrals con¬ 
verge simultaneously. For instance, the energy £2 is finite for n > 1 
if fll = 0 or for n < — 1 when 03 = 0. The energy £3 is given by the 
formula 

E 3 = l)(2ai -|-a2+ 2a3) + (ye^ - l)^^] 

where the coefficient 04 reads 


fl'^ = 32{Pe^ + ye'^-2)n^K^cr- P[ cf. 


«4 E 
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(« 1 , 

« 2 ) 

£ 

£1 

£2 

£3 

£4 

£5 

£6 

El 

(3, 

1 ) 

224 

18.9 

70.0 

86.4 

22.6 

0.23 

2.93 

22.7 

( 2 ,- 

- 1 ) 

226 

18.9 

71.7 

86.5 

22.6 

0.22 

2.97 

22.8 

(4, 

2 ) 

321 

25.4 

84.5 

127 

34.3 

0.60 

15.2 

34.5 

(4, 

1 ) 

276 

25.3 

81.4 

112 

31.4 

0.59 

-6.46 

31.8 

(3,^ 

- 1 ) 

277 

25.3 

82.1 

112 

31.3 

0.59 

-6.35 

31.9 

( 2 ,- 

- 2 ) 

328 

25.3 

90.4 

127 

34.2 

0.59 

15.5 

34.6 

(5, 

2 ) 

378 

31.5 

93.4 

158 

42.7 

0.28 

9.40 

42.5 

(5, 

1 ) 

325 

31.8 

89.8 

126 

41.8 

1.53 

-8.99 

42.9 

(4,- 

- 1 ) 

324 

31.8 

89.2 

126 

41.8 

1.51 

-8.76 

42.9 

(3,- 

- 2 ) 

379 

31.5 

94.8 

158 

42.5 

0.28 

9.78 

42.7 


Table I: The energy and the components (in unit of AM^) with 
the potential (l46ll,(l47ll of {a,b) = (0,2). The parameters are 
(i3e2,ye2,^2,fei,A:2) = (2.0,2.0,1.0,1.0,2.0). 


The term £4 takes the form 


„4 


£4 = 


-M^(/3e^ + ye^- 2 )c' 
i \n\ 


'N-K+\ 


16%^ 


■M^{jie^ + re -2 )cT^ 


K 

N-K+\ 


K 


- l-i 

7o -''(iTtr 

« 2 -l 


^”(r) 


This contribution to the energy does not appear in the models without 
potential where /3e^ + ye^ — 2 = 0. The terms do not con¬ 

tribute to a total energy since they vanish for a holomorphic solution. 
The last contribution which comes from the potential term reads 


£7 = 


rdppv 

Jo 




[^] 


£2(iV-£ + l)2 Jo 


" r^^TT 

Jo 1 




where the function Q is taken for \uj'\^ = 1. Considering that p? 
is given by ( 128b one gets that £7 = £ 4 . It follows that for an exact 
solution the potential term and the quartic term proportional to jSe^ + 
ye^ — 2 contribute equally to the total energy. 


IV. THE NUMERICAL STUDY EOR THE 
NONHOLOMORPHIC VORTICES 


In present section we study the problem of solutions of the ex¬ 
tended CP^ Skyrme-Faddeev model without reduction and with 
presence of the potential. We shall propose the form of the poten¬ 
tial such that one can compute solutions for arbitrary set of integers 
{n\,n 2 , ■ ■ ■ pin) which appear in the ansatz ([8l(. Such solutions are 
non-holomorphic ones and therefore they can be obtained as the re¬ 
sult of numerical integration. For the numerical study, it is more 

convenient to use a new radial coordinate y, defined by p = yj 
Accordingly we adopt profile functions gi, instead of using /,. The 
ansatz is then 


^ Si{y) j(m (p+kj y/'fw)) 

gi{y) 


(33) 


where vr(w) is a real function of the light-cone coordinate w = z ± f. 
The factor I/a/V is introduced in order that the solution naturally 
reduces that of CP^ when all integers n/ are equal. It is worth men¬ 
tioning that in this section we are not interested in reduction itself, 
however, for the case of reduction one can test the numerical solu¬ 
tion comparing it with analytical one. 

The equation of motion ([IB can he written as 


■-f I 


c'l , i- 2 y 


■-f 


): 


8i 




g,(l-g/) VCi y(l-y)/g;(l-gi) 2 (g,(l-g ,))2 

1 /C( 




Vy(l-y)^ VCi 2 Ciy(l-y) 


1 \ 1 1 

l-y))”'’‘' 2 v(l- 


C3 2 
n: 


y(l-y)3Ci 


0 

P^ 1 




Si 


1 


C 3 


^iTi^gjJ8i{i-8i) y( 1 -y)^ Cl g 

0 


8 y3(i^3,)Ci 
where 0 = 8k ^ symbols Q read 


5V,=0, i=l,---,A (34) 


Cl = Cpp — —1 4 - ( 2 ye^ +/3e^)6(y) -t- {jje^ — 1) 

Cl = iCp(f = —3^(y) 

C'i, = C(p(p 

= -!^ + (iSR _ i)i^e(y) + (2ye2 + i3R)®(y) (35) 

y y 

and dVi stands for contribution from the potential i.e. 


, SV 


SV 


5V = -Ag 2 —-^g,(l-g,) —. 

^Si ^Sk 


(36) 


There is a freedom in choice of the form of the potential V (g,), 
however, one has to take care about its asymptotic behavior. A stan¬ 
dard discussion should be based on the vacuum structure of the field. 
As a simple example, we start with the C£* (0(3)) case. The 0(3) 
C7 model is usually defined as a vectorial triplet n= {111,112,112) with 
the constraint n-n= 1. The well-known potential named “old-baby” 
type, i.e.potential with one vacuum, is of the form 


V (i?) = (1 — iioo • n) 


(37) 


where ra„o is a vacuum value of the field H at spatial infinity. If we 
choose the value Hoo = ( 0 , 0 , 1 ) the potential becomes V = (1 — 173 ). 
Performing stereographic projection on a plane we parametrize 
the model by a complex scalar field u related to the triplet n by 


n =- T—pr (u + u*,—i{u~u*),\u\^ — l') 

l + |«p \ / 

and rewrite the potential in terms of a complex field u 

2 


V(m)= 1- 


1 + IkP 


1 + m 


2 • 


(38) 


(39) 


One can expect that similar argument might work for the CP^. In or¬ 
der to check this hypothesis let us consider two fields (mj , M 2 ) whose 
behavior at the infinity is the same as for corresponding holomorphic 
solutions characterized by two integers (ni,n 2 )- When ni,n 2 > 0, 
the fields behaves as |mi|,|m 2 | ^ °° for p —> Then one can try a 
generalization \u\^ -P -f |m 2 |^ resulting in the potential 


V{ui) 


|MlP + |M2p- 1 \ _ 2 

i + |Mip + |M2py 


1 

Uj = 

yfN 


l + |MlP + |M 2 p' 


(40) 
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Figure 1: The plot of the CP^ profiles gi (the solid line) and g 2 (the dotted line) for the potential of {a,b) = (0,2). The parameters 

are = (2.0,2.0,1.0). 


There is a serious problem which such a generalization since the 
model with the potential l l40b has numerical solutions only for equal 
values of integers tii = n2- 

A better approach to the problem is based on the observation that 
the potential in the CP^ case can be expressed in terms of the SU (2) 
valued field U :=z-n which allows to write l l43b as 

V{U) = ^Tr{l~UlU). 


This formula can be easily verified using the identity (?-nc«,)(?-M) = 
Hoo • n + i7 • (Hoc X ii). We examine the construction of the potential 
for the A > 1 of the CP^ case in the same way. A parametrization 
of the model, which includes the case of SU(2)/U{1) = CF* target 
space, is performed in variable X instead of the St/(2)-valued field 
U. The principal variable is a function of one complex field u and it 


(41) 
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Figure 2: The plot of the CP^ density of total energies ji" (in unit of 4M^) of the {k\,k 2 ) = (1.0,1.0) (the solid line), the {k\,k 2 ) = (1.0,2.0) 
(the dotted line) and (fci, kx) = (2.0,1.0) (the dot-dashed line) for the potential ( I46b .( l47b of {a, b) = (0,2). The parameters are (/3e^, ye^, /t^) = 
( 2 . 0 , 2 . 0 , 1 . 0 ). 


reads 

XCP' - ^ \ 

1 + |MP \ liu* l-|Mp/ 

or in terms of components of the unit iso-vector n 

= —tr'I + m' Ti — in^Z2, 


which is clearly different form U. Inverse of the principal variable 
goes to (X^^')^* = diag(—1, — 1). It follows that the following 
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potentials of which the solutions share the asymptotic behavior with 
the holomorphic counterpart, i.e. ~ (p”‘ ,p"2) for the combinatation 


1. The case: n\ >112 >0 

By assuming that the solution and its holomorphic counter¬ 
part have the same asymptotic behaviour at the spatial infinity 
one gets that inverse of the principal variable X goes to := 

diag(—1,1, —1) as p —t 00. It follows that generalization of the for¬ 
mula I l43l l from = 1 to = 2, gives the following expression for 
the potential 


V(u;) = Tr(l-Xco 


^X)=4 


1 + 

l + |MlP + |i<2p 


(45) 


Figure 3: The total energy surface in the parameter space (j8e^,7e^) 
potential of {a,b) = (0,2). The remaining parameters are 

(A2,fci,fc2) = (i.o,i.o,i.o). 


Note that for p —t 0 inverse of the principal variable goes to 2 lq * : = 
diag(l, 1,1), then the expression Tr(l — can be included as 

the “new-baby” potential which has two vacua I14ll . Finally, the fol¬ 
lowing expression can be considered as a general form of the poten¬ 
tial 


expression 

1/ = [Tr(l-AQ"'A)]“[Tr(l-A,;;'A)]* 


V{u) = ^Tr{l-{xCP')-^X^P') = ^Tr{l+xCP') 

(|ui|2 + |u2 |2)«(1 + |u2 |Y 


(1 + Ulp-f M2p)“+* 



(gl+g2 - 2gig2)‘'gt(l+g2)'' 
(gl+g2)“+'' 

(46) 

reproduces the potential ||39I. The last result constitute an important 
clue how to choose potentials V (A) for A > 1. 

where the integers a, b satisfy a ^ 0,b > 0. 



A. The CP^ solution 


2. The case: n\ > 0> 112 


First we give a definition of functions which still were not ad¬ 
dressed. For the CP^, the functions 9,co,l^ in EQ.l l35t have the form 


eiy) 



(1 -t) 


' (l+g2)g2gf 

. 2 gi(l-gi) 


(l+gl)glg? 
2g2(l -g2) 


-g'lg'l 


0 }{y) 



«1 gl (1 - gl )g2 + «2glg2 (1 - g2) 


+ 2("1 -«2)^gl(l-gl)g2(l-g2) 


Assuming that for U2 < 0 the field U 2 behaves at zero as its 
holomorphic counterpart i.e. ~ p"^ one gets that it tends to di¬ 
verge as p —!■ 0. Then inverse of the principal variable X goes to 
:= diag(l, —1, —1) as p —S’ 0. The general form of the potential 
takes the form 


^ (1 + | ui | 2)«(1 + | u 2 |)^ 
(l + |uip + |u2p)“+* 
^ gtg2(l+gl)°(l+g2)" 

(gl +g2)“+'' 

where the integers satisfy a ^0,b > 0. 




nigh'l +«2glg2 


+ 2(^1 -"2){g2(l-g2)gi -gl(l-gl)g2} 


(44) 


3. The case: n\,n 2 <0, |«i| > |n2| 

The asymptotic values of inverse of the principal variable 
are given by constant matrices = diag( 1,1,1) and X^^ = 
diag (— 1,1, — 1). Then the potential is 


where © = gi -l-g2- 

Generally speaking, a potential can be deduced from the asymp¬ 
totic structure of solutions of the model. Moreover, the potential has 
to have the form such that the model has solutions for all qualitatively 
different combinations of the integers (ui,«2)- In the following sub¬ 
sections, we give an explicit form of the potentials for basic combina¬ 
tions of {ill , 112 ). The most crucial point is that we shall explore such 


^ (l + |u2p)‘'(|niP + |u2p)'' 

(l + |i<lp-|-|M2p)“+* 

gl(l+g2)‘'(gl+g2-2gig2)'’ 

(gl+g2)“+'' 

where the integers satisfy a^0,b > 0. 


(48) 
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Figure 4: The profiles (gi,g 2 ), and the energy density of the CP^ in the potential l l46t of (a^h) = (0,2) with several strength For the 
remaining parameters we choose (j3e^,ye^,^i,fc2) = ( 2 . 0 , 2 . 0 , 1 . 0 , 2 . 0 ). 


There still be a freedom for choice of the parameters {a,b). Here 
we consider the simplest case, i.e., (a, b) = ( 0 , 2 ) for which the forms 
l l46l l and l l47t become identical. (Note that from the asymptotic anal¬ 
ysis there are no solution for b= \.) From i36t . the contributions of 
the potential term 5V/ can easily be estimated as 




. g?g2(l+g2) 
(gl +g2)^ 


5V2 = 0. 


(49) 


Expansions 

We examine the asymptotic behavior of the solutions expanding 
the equations with l l49b . the asymptotic behavior at the spatial infinity 
p —>■ oo (y = 0 ) is given by the series expansion 

gi(y) = ciy"‘+0(y"‘+') 

g2(T) = c2/^+0(/^+‘) (50) 

for ,«2 >0 and 

gi(y)=ciy”>+0(y”' + ') 

g2(y) = l+C2y"^+0(/^+‘) (51) 


has logarithmic divergence at y = 0. For the cases such as (n i, n 2 ) = 
(3,1) the density is 


(3,1) _ 8ci(*i -fc2)^„2 , 


> = 


C2 


-y^ + Oif). 


(54) 


One can easily see that it leads to finite energies per unit of length. 
In the following part, we will present results only for non-divergent 
cases. 

The numerical analysis is performed by a standard relaxation tech¬ 
nique of which a typical mesh size is chosen as Al„,esh = 1000 , which 
supports the good convergence property of the solution. In Table 
m we summarizes values of the energy and the components. The 
El is the topological term, i.e. its value is 2 ;r«i (for «2 > 0 ) or 
2;r(«i + |u 2 |) (for >^2 < 0). Our results are qualitatively good and 
the uncertainty is less than 1 percent. Again note that the Derrick’s 
scaling argument for two spatial dimensions implies that the energy 
per unit length from the quartic terms and the potential terms should 
balance,!.e., £4 +£5 = £ 7 . Then, the maximal vlaue of the uncer¬ 
tainty of our numerical results is ~ 1 percent. Fig[T]plots the several 
profile functions. Now we define the hamiltonian densities in 
terms of the energy per unit length 


for «i > 0,«2 < 0. The c, is arbitrary constants (“shooting param¬ 
eters”) and then all the higher order coefficients can be written by 
Cj. It has been also checked for the present potential that expanded 
solution has good asymptotic behavior at the origin p —>■ 0 (y —>■ 1 ). 

The fact that our numerical solutions and holomorphic solutions 
have the same leading asymptotic behavior means that they share the 
problem of convergence of the energy contributions. The analysis 
performed in H shows that not all combinations of integers (ni ,^ 2 ) 
leads to finite energy per unit of length. The most troublesome term 
is For this reason we shall study its asymptotic behavior. For 
instance, the term has the leading expansion term around y ~ 0 
(for the case of V^(w) = w) 

(2.1)^ 8ci(fcl-k2)^ y^ 2) (52) 

C 2 

which causes the divergence of the integral unless ki = k 2 because 
the integral 

J= 2k ^’^ 2 ^’^'’ (53) 


Ek = 2 kM^ [ pdpMi = Aiyf P (55) 

Jo Jo 2y^ 

In Figl^we present the JiP as a function of radial coordinate p for 
some values of (^ 1 ,^ 2 )- Fig [3 shows the total energy surface in the 
model parameter space j3e^ and ye^. 

When pe^ + ye^ = 2, the holomorphic solutions 19} are scale in¬ 
variant, then the coefficients are freely chosen. The fourth order 
terms proportional to i.e. pe^ + ye^ 7 ^ 2 , breaks the scale invari¬ 
ance of the model what leads to the fixing of the coefficients at some 
values. It turns out that the presence of only such fourth order terms 
does not lead to numerically stable solutions. In order to find so¬ 
lutions, we introduce the potential which fixes the solution corre¬ 
sponding to the highest integer n^. As a consequence, the Derrick’s 
theorem is satisfied because the coefficients of all remaining compo¬ 
nents are properly determined. It is worth to examine the behavior of 
the size of solutions as a function of the strength of the potential p^. 
In Fig|4l we show the plot; one can see that when increasing the 
solutions became better localized around the center. 
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Figure 5: The profiles (gi,g 2 :.? 3 )i ™d the energy, their components of the CP^ in the potential ( 159b of {a,b) = (0,2). For the remaining 
parameters we choose {j5e^ ,Ye^, p.^ ,ki,k 2 ,k^) = ( 2 . 0 , 2 . 0 , 1 . 0 , 1 . 0 , 2 . 0 , 1 . 0 ). 


B. The higher N solutions 


V. SUMMARY 


A generalization of the presented approach for the A > 2 is almost 
straightforward. Now we consider the case of all integers are posi¬ 
tive, i.e. > 0,t = I,-" ,N. If Ml has the highest positive integer, 
inverse of the principal variable X goes to Aoo = (—I,!,-- - ,1,-1), 
then the general form of the potential is then 


V{ui) 


iLl, I 


(1+I7=2l 


\ 2 \b 




(56) 


with fl ^ 0 , > 0 . If the two highset positive integers are equal = 

n 2 , the potential naturally reduces to the We examine the 

case of the CP^. Inverse of the principal variable X goes to 


/ 0 -10 0 \ 
^_l^ -10 0 0 

0 0 10 
Vo 00-1/ 

which results in form of the potential 

V{ui) = Tr(l-Aj‘A) 

^ 2-F|mi -M2p + 2|M3|^ 

l + |Mip-|-|M2p + |M3p' 


(57) 


(58) 


When we put u\ =U 2 = uj'Jl., the result coincides exactly with the 
result obtained for CP^ E5}. Here we examine the case of which 
n\ > n 2 > 112 > 0 and the expressions SVi have the following form 


s,, ,g?g2g3(g2+g3+g2g3) 

OVl = —O—. - r~- 

(8283 + 8381 +8l82r 

SV2 = SVs = 0 . (59) 

A typical result of the CP^ is shown in Fig[3 The value of the topo¬ 
logical term is Ei ~ 25.3 and the combination of energies per unit 
length concerned with the Derrick’s becomes E 4 + ~ ~ 0.12. 


In the present paper we were focusing on the problem of solutions 
of the extended CP^ Skyrme-Faddeev model for wide range of cou¬ 
pling constants i.e. when + ye^ 7 ^ 2. The results of our analytical 
and numerical studies indicates that such solutions do exist, however, 
they require a presence of the potential term. In the first part we have 
considered some special cases when the reduction from the CP^ to 
the Cf' happens. A particular choice of the potentials in such cases 
enabled us to obtain the holomorphic solutions in the sector of cou¬ 
pling constants 7 ^ 2. There is still an open problem if such 

solutions do exist for non-reduced CP^ case. 

For a general case a freedom of choice of the potential is much 
larger. As a consequence our potentials consist of one vacuum type 
(“old-baby” potential) as well as of two vacua type (“new-baby” po¬ 
tential). Both potentials break the scale invariance of the solutions 
and then they satisfy the Derrick’s theorem. The numerical solutions 
presented in the paper does not satisfy the zero curvature condition 
and therefore they are not holomorphic functions, however, they have 
the same asymptotic behavior as its holomorphic counterparts char¬ 
acterized hy nj^. 

In this paper, we have numerically examined only the solutions 
with the rotational symmetry in the case of the old-baby type poten¬ 
tial. The new-baby potential should have the similar solutions, too. 
It is a well-known fact that the old-baby potential tends to split solu¬ 
tions with topological charge B into B independent fractions, while 
the new-baby type does not It would be interesting to inves¬ 
tigate if such behavior do really manifest for non-central vortex so¬ 
lutions in the model with our potentials. Such a study is important 
since it can serve as a test for validity of our proposal for that poten¬ 
tial. The analysis of this subject is in progress and the results will be 
reported in near future. 
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